Effect of electron-phonon interaction on the impurity binding energy in a quantum 

wire 
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The effect of electron-optical phonon interaction on the hydrogenic impurity binding energy in a 
cylindrical quantum wire is studied. By using Landau and Pekar variational method, the hamilto- 
nian is separated into two parts which contain phonon variable and electron variable respectively. 
A perturbative-variational technique is then employed to construct the trial wavefunction for the 
electron part. The effect of confined electron-optical phonon interaction on the binding energies of 
the ground state and an excited state are calculated as a function of wire radius. Both the electron- 
bulk optical phonon and electron-surface optical phonon coupling are considered. It is found that 
the energy corrections of the polaron effects on the impurity binding energies increase rapidliy as 
the wire radius is shrunk, and the bulk type optical phonon plays the dominant role for the polaron 
effects. 

PACS: 71.38+i;73.20.Dx;63.20.Kr 



I. INTRODUCTION 



During the past decades the development of the epi- 
taxial crystal growth techniques such as molecular beam 
epitaxy and metal-organic chemical vapor deposition has 
made the growth of the quasi-two-dimensional (quantum 
well) or quasi-one-dimensional (quantum wire) Q |^ |^ Q 
systems with controllable well thickness or wire radius 
became possible. These quantum structures have 
been applied to many semiconductor devices, such as 
high-electron- mobility transistors. Recent progresses in 
growth and fabrication techniques have been able to 

o 

fabricate the quantum wires with radii less than 100 A. 
Theoretically, the electronic properties of a hydrogenic 
impurity in the quantum well[a,S0la| and the quantum 
wire jHHIil m m m have been studied by many 
authors. The impurity binding energies of a quantum 
wire with infinite or finite potential barri er |g and with 
different shapes of the cross-secti on|ip[ have been 
discussed. The effect of locationjlfl [HI of impurities 
with respect to the wire axis was also studied previously. 
The emission line for quantum wires was observed Tsj 
to be two to three times broader than that of quantum 
wells and with 6-10 meV higher binding energy. It is 
expected that the same properties in quantum wells were 
further improved by the reduction of dimensionality to 
quasi-one-dimensional quantum wires. 



The physics of impurity states in quantum wire is 
very interesting because specific properties can be easily 
achieved by varying the wire radius. An electron bound 
to an impurity on the axis of the quantum wire behaves 
Hke a bounded three-dimensional electron when the 
boundary is far away. However, as the wire radius is 



reduced, the electron confinement due to the potential 
barrier becomes very important. Especially in the 
quantum wire with infinitely high potential wall, the 
total energy of the electron may change from negative to 
positive at a certain radius and finally diverges to infinity 
as the radius approaches zero. Furthermore, it is well 
known that the reduction of dimensionality increases 
the effective strength of the Coulomb interaction. The 
binding energy Eb of the ground state of a hydrogenic 

impurity in N-dimension is given by E},— 



R*. 
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where i?* — jpj^ is the effective Rydberg. Hence the 
dramatic change in the binding energy may serve as a 
clear signal for variation in the effective dimension of 
the quantum wire. 



It is known an electron weakly bound to a hydro- 
gen impurity in a polar semiconductor will interact with 
the phonons of the host semiconductor. In the past 
decade, many authors have studied the polaron effect 
on the binding ener gy o f impurity or exciton in quan- 
tum wellEl 111 Ellli m In il El H. Recently, 
the electron-phonon effect on the binding energy of the 
donor impurity in a quantum wire with rectangular cross- 
section was reported[2l It was found the po- 
laron effect on the binding energy becomes sizeable as 
the electron gets more deeply bound. The polaron shifts 
in donor energy levels are found to be of the order of 
10% in a weakly polar system. In studying the polaron 
effect on the impurity binding energy, most of the previ- 
ous works considered the interaction of the electron and 
bulk optical(BO) phonon only. However, in ionic crys- 
tal, the motion of an electron near the surface may be 
affected very much by the surface longitudinal optical 
(SO) phonon [2^. An electron may be trapped at the 
surface by the electron-SO phonon interaction. Besides, 
the electron phonon interaction Hamiltonian in the previ- 
ous works was valid only for the bulk. Therefore, we will 
choose the Hamiltonian derived by Li and Chen [23], who 
considered the confined phonon modes in the cylindrical 
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quantum dot. 

Most of the previous approaches concentrating on 
the polaron effect on the ground state of an impurity 
in a quasi-one-dimensional wire employ the variational 
method or perturbation method. Since the construction 
of variational trial wave functions bases entirely on phys- 
ical intuition, and the estimation of the accuracy of the 
result obtained from variational approach is very diffi- 
cult. Furthermore, the perturbation method is only a 
good access to those systems with very small perturba- 
tion in most cases. Therefore, it would be most desir- 
able to have an alternative approach which is not only 
simple but also efficient to the quantum wire problem. 
In this work, we employ a simple approximation treat- 
ment which combines the spirit of both variational prin- 
ciple and perturbational approach to study the effect of 
electron-phonon interactions on the ground state binding 
energy of a hydrogenic impurity located inside a quantum 
wire. 



= E-^o(^p)[ E Vnicos{^z){a^i+al) 



1=1,3,. 



2D 



Itt 



+ Yl ^»'Sin( — z)(a„; -|-a,\;)] 



1=2,4,. 



with 



(5) 



V 



ixJdfJHxJ + ii^nDrJKxn) 



2 72/ 



1 

(6) 

where Jm is the mth-order Bessel function , Xn the 
nth-root of Jq , and V = 2Tr(PD{D — > oo) is the crystal 
volume. Hsp is the surface optical phonon (SO) phonon 
Hamiltonian which can be expressed as: 



II. THEORY 



(7) 



Consider now a hydrogenic impurity located on the 
axis of a rigid wall cylindrical quantum wire with a radius 
d. The Hamiltonian of the impurity electron interacting 
with the phonon can be expressed as: 



where hojsp is the surface optical (SO) phonon energy, 
Bl^{Bn) is the creation (annihilation) operator for SO 
phonon. i?e-sp is the interaction between electron and 
SO phonon: 



H — He + Hb + He-b + Hsp + He-sp, (1) 

here, H^. is the electronic part of the Hamiltonian 



He 
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+ V{p) (2) 



2/i \dx^ dy^ J 2^dz^ er 
V{p) is the confining potential which is assumed as: 

for p < d, 



V{p) = 



00 for p > d, 



(3) 
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and e and fi are the dielectric constant of the well and the 
effective mass of the electron. Recently, Li and Chenp9l| 
has derived the confined the longitudinal-optical phonon 
and surface phonon modes of a free-standing cylindrical 
quantum dot of radius d and height 2D. We will fol- 
low their Hamiltonian and let D approach infinity, such 
that the dot system can become a quantum wire. There- 
fore, Hi, is the bulk phonon Hamiltonian which can be 
expressed as: 



n.l 



(4) 



where Tiujlo is the dispersionless bulk optical (BO) 
phonon energy, a|jj(a„;) is the creation (annihilation) op- 
erator for BO phonon. He-b is the interaction between 
the electron and BO phonon which can be expressed as: 
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where kn — , and S — nd^ . Im and km are , respec- 
tively , the mth-order modified Bessel function of the first 
and second kind. 

Following Landau and Pekar's variational 
approach (SOj. the trial wavefunction can be written 
as: 



^' > = Hr)UbUs I >, 



(12) 



where $(r) depends only on the electron coordinate, and 
I > is the phonon vacuum state defined by 6g | >= 0, 
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>= 0, and U is a unitary transformation given by: 



Ub = exp ^{ali fni - ani f* 



(13) 



(14) 



Where /„; and gn are the variational function and the 
unitary operators Ub and Us transform the bulk phonon 
and surface phonon operators as follows: 



UlunlUb = ani + fnl, 

UtBiUs = Bi 
UlB^Us = B„ 



gn- 



(15) 
(16) 
(17) 
(18) 



The parameters /g, /*, gq, can be obtained by mini- 
mizing the <\ H \> with respect to the parameters 
/*;, gm g*. Then the expectation value < H > turns out 
to be 



of the Hamiltonian He contains two terms, i.e. Hqi and 
Hq2 , where Hqi represents the one dimensional harmonic 
oscillator, and H02 represents a two dimensional hydro- 
gen atom located inside a quantum disk|l4j. Both can 
be solved exactly. For illustration, the ground state en- 
ergy and wavefunction of the unperturbed part can be 
expressed as : 

4"i)(a,/3) = 4"i)(/3)+£;f)(a) (24) 
^<f\r,a,P) = )(^;/3)¥.f )(x,2/;a) (25) 

respectively, where (p^^^\z; P) is the ground state wave- 
fimction of the ID harmonic oscillator, and (p^^^\x, y; a) 
is the ground state wavefunction of the 2D hydrogen 
atom located at the center of an infinite circular well. 
The ground state eigenvalue and eigenfunction of the ID 
harmonic oscillator can be expressed as: 
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(27) 
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The ground state eigenvalue and eigenfunction of the 2D 
hydrogenic impurity located at the center of an infinite 
circular well can be obtained asjQ: 
(1) For E <0, 



< *(r) I /o(fc„p)cos( — z) I $(r) > 



The axis of the wire is assumed to be along the z di- 
rection. To solve the electronic part, one can employ 
the perturbative-variational approach as we did in the 

above subsection. Two variational parameters a and /3 

2 

are introduced by adding and subtracting two terms 



and 
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into the original Hamiltonian H^ and then 



regroup H^ into three groups: 



He^Hoi{P)+Ho2{a) + H'{a,(3) 



where 
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(21) 
V{p) (22) 
(23) 



In the above equations, H'{a,/3) is treated as a pertur- 
bation, and a and P are treated as variational parame- 
ters which can be determined by requiring the perturba- 
tion term to be as small as possible. Decomposing H^ 
into two terms Hqi and Hq2 is equivalent to dividing the 
space into a two-dimensional (in xy plane) and a one- 
dimensional (in z-axis) subspace. The unperturbed part 



ipf^\x,y;a) = Nie~^^["'^^^<^{\ m |+l/2-Ai,2| m 

(28) 

where = a^p, ai = Ai = ff^, $(a, b, x) is the 

confluent hypergeometric function, and A'^i is the normal- 
ization constant. 
(2) For E>Q, 



<^f)(x,y;«) = iV2er<i>„ 1(7/2,^2) 



(29) 



-J^^: ^m-l/2('72''C2) 



where ^2 = "2^, ^2 = —jf^, V2 

is the irregular Coulomb wave function, and N2 is the 
normalization constant. 

(3) The turning point for energy changing from E > to 
-E < in the quantum circle system may be determined 

by setting 



8pe'^ 



1/2 



for m = 0, (30) 



and 



d-'/'j2 



8^2 X 1/2 



,1/2 



for m = 1, 



(31) 



The requirement of the continuity of the wavefunctions 
and its first derivative at boundary yields: 
(1) For E <0, 

(^(1 m I +i-Ai,2| m I +l,a,d) = (32) 
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(2) For £ > 

^m-l/2iV2^a2d) = 

The eigenvalues are then given as: 



-2fi% ,forE<0, 
#l4 '/«'^^>0' 



(33) 



(34) 



The first order energy correction can thus obtained as: 
A£;W(a,/3) = <^f\r;a,(3)\H'ia,p)\^f\r;a,p)> 



-<<i)(^;/3) 



ep 

2jLl 



^f)(x,y;a)> 



-<$(°)(r;a,/3)|-|$(")(r;a,/3)> 



The second term of the above equation can be integrated 
analytically and the result is: 




FIG. 1: The ground-state energy(solid line) and the hinging 
energy (dotted line) of a hydrogonic impurity located at the 
axis of a cylindrical wire as a function of the radius of the 
wire. (Ry* and a* are effective rydberg and the effective Bohr 
radius.) 



2/i 



(/pf)(^;/3)>^f)(x,y;a) >= 



fj%^ and without the impurity, i.e. 



(35) 

Then the total energy up to the first order perturbation 
correction can then be obtained as: 

E,{a, p) = Ef^^ (/3) + 402) + ^^(1) (c,, /J), (36) 

The variational parameters are then chosen by requiring 
the total energy Eg{a, (3) to be minimized with respect to 
the variation of a and (3. This is equivalent to requiring: 



dE 
da 



dE 
'80 



= 0. 



(37) 



(38) 



For the excited states, the eigenvalues and eigenf unctions 
can be treated in the same way. 



III. RESULTS AND DISCUSSIONS 

We have calculated the effect of the confined the 
longitudinal-optical phonon and surface phonon interac- 
tions on the hydrogcnic impurity located in a quantum 
wire. And the well potential is considered as infinite. 
Figure 1 shows the ground state energy as a function of 
the wire radius. The binding energy £{, of the hydrogenic 
impurity is defined as the energy difference between the 
ground state energy of the cylindrical wire system with 



Eb = Eq — Eg 



(39) 



where Eq is the ground state energy of the quantum wire 
system without the impurity, while Eg is the ground state 
energy of the quantum wire system with the impurity lo- 
cated on the axis of the cylindrical wire. One can see 
from Fig.l that the energy of the Is state becomes nega- 
tive when the wire radius is larger than 1.65a*. It means 
that the confining energy is larger than the Coulomb en- 
ergy as the wire radius is smaller than 1.65a*. And one 
can also note that as the radius of the quantum wire 
is decreased, the ground state energy increases. As the 
wire radius d becomes smaller, the electron is pushed to- 
ward the axis of the cylindrical wire. This makes the 
electron get close to the nucleus. As the electron gets 
close to the nucleus, both the ground state energy and 
the binding energy increase rapidly. This is because the 
Coulomb potential, which varies with ~ i {d is the wire 
radius), becomes more negative, while the kinetic energy 
of the electron, which varies with (by the uncer- 

tainty relation), increases more rapidly. As a result, the 
ground state energy is increased as the electron gets close 
to the nucleus. The binding energy defined in Eq. (39) 
is effectively the negative sign of the of the Coulomb in- 
teraction energy between the electron and the nucleus, 
i.e. ~ ^, therefore, the binding energy of the electron is 
also increased as the electron gets near to the nucleus.. 
As a result, the groimd state energy is increased as the 
electron gets close to the nucleus. The binding energy 
defined in Eq. (39) is effectively the negative sign of the 
of the Coulomb interaction energy between the electron 
and the nucleus, i.e. ~ ^, therefore, the binding energy of 
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the electron is also increased as the electron gets near to 
the nucleus. Our results show that for small wire radius, 
the binding energies are in good agreement with previous 
results [TlLll4j . As the radius becomes very large, our re- 
sult ajpproaches the correct limit IR* while the previous 
work can only yield a value of 0.22i?*. The large 
discrepancy of the previous work may be due to the ar- 
tificial dividing of the variational trial wavefunction into 
a one-dimensional hydrogen atom and a two-dimensional 
hydrogen atom and thus forces the creation of an addi- 
tional node of the wavefunction at z=0. In this work, the 
trial wavefunction is adopted to be the form of ID har- 
monic oscillator wavefunction instead of one dimensional 
hydrogen atom. This prevents our wavefunction from in- 
troducing any additional node at z=0. Figure 2 presents 
the 2s excited state binding energies as the functions of 
wire radius. One can note from the figure that as the wire 
radius increases, the binding energy approaches 0.25i?* 
which gives correctly the limiting value of 3D hydrogen 
atom. 



— t:onfiiied BO-phoiion 
confined SO-phonon 
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wire nidiusla*) 



FIG. 3: The energies of the wire radius modified by the con- 
fined BO and SO phonon. The solid line stands for the BO 
phonon effect, and the dashed line for the SO phonon effect. 
(Ry* and a* are effective rydberg and the effective Bohr ra- 
dius.) 
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FIG. 2: The binding energy of the 2s excited state as a func- 
tion of the wire radius. (Ry* and a* are effective rydberg and 
the effective Bohr radius.) 



face area of a quantum wire (or quantum dot) decreases 
with the radius. Thus the number of vibration modes of 
confined SO phonon becomes fewer. 
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Figure 3 presents the confined BO phonon and SO 
phonon effect as a function of wire radius. With increas- 
ing the wire radius, the magnitude of the confined BO 
phonon effect decreases from large value and then ap- 
proaches to the bulk value. When the wire radius is less 
than 1.5 a*, the polaron effect increases rapidly. One 
might think as the radius becomes very small, the con- 
fined BO phonon effect should approach zero, like the 
case in quantum well [sj. In fact, similar results were 
obtained by Oshiro in a spherical quantum dot [3.^. They 
found the polaron energy shift is enhanced as the dot ra- 
dius becomes small. This is due to the fact that the 
electron becomes complete localized (E;, approaches in- 
finity) in small wire (or dot) radius while the binding 
energy approaches 4i?* in small well width. In the case 
of quantum well, the confined SO phonon effect plays the 
dominant role for small well width . But in quantum 
wire, the confined SO phonon is less important, just like 
that in quantum dot system "s^l . This is because the sur- 



wire radius(a*) 

FIG. 4: The binding energy with/ without phonon effect. The 
dotted line stands for the binding energy without the phonon 
effect. The dashed line stands for that only BO phonon effect 
on the binding energy, and the solid line for both the BO and 
SO phonons effects on the binding energy. (Ry* and a* are 
effective rydberg and the effective Bohr radius.) 

In Fig. 4, three curves are presented. The dotted curve 
represents the binding energy of the impurity without 
considering the interactions between the electron and 
phonon. The dashed curve represents the binding energy 
of the impurity with only confined BO phonon effect be- 
ing taken into account. While the solid curve is the bind- 
ing energy of the impurity including both confined BO 
phonon and SO phonon effects in the calculation. Com- 
paring to the impurity binding energy, the confined SO 
phonon is negligible in quantum wire. We then conclude 
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that because of the similarity in geometry, the behavior 
of the polaron effect on the quantum wire system is hke 
that on the quantum dot system. 

IV. CONCLUSION 

In this work, analytical solutions for the effects of the 
electron-phonon interaction on the binding energies of an 
impurity located inside a quantum wire are obtained by 
a simple but efficient perturbation-variation method. As 
the radius becomes very large, the correct limiting value 



can be obtained. Wc have also discussed both the con- 
fined BO and SO phonon effects. We found the confined 
BO phonon effect is prominently for a quantum wire with 
small radius. We also found that the energy corrections 
of the polaron effects on the impurity binding energies 
increase rapidly when the wire radius is less than 1.5 a*. 
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